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Abstract. Starting for the Stillinger and Weber expression for the free 
energy of supercooled hquids, we extend the free energy to the case in 
which two time scales separate and the system is in quasi-equilibrium. The 
concept of an effective temperature, different from the kinetic temperature, 
is naturally introduced. An example of hypotetical quasi-equilibrium phase 
diagram is presented for the case of SPC/E water. 



1. Introduction 

The study of the properties of undercooled liquids is a topic of current 
research. Under the name of undercooled liquids we will indicate liquids 
whose relaxations start to show a separation between different timescales. 
The idea of considering the dynamics of undercooled liquids as similar to 
the vibrations of a disordered solid on short timescales, while attributing 
the slow dynamics to rare change in the global structure (i.e. the atomic 
position of the disordered solid assumed as the reference for short time 
vibrations) is old [1]; in particular, the fruitful reformulation of Goldstein [2] 
of these ideas in terms an energy landscape has led to the introduction of a 
formalism [3] that is particularly suitable for the study of liquids via Monte 
Carlo or Molecular Dynamics simulations [4]. 

In Sec. 2 we describe such formalism, showing that it is particularly suit- 
able to investigate the connections between the statics and the dynamics of 
equilibrium systems. In Sec. 3 we show that such formalism can be extended 
to systems in quasi-equilibrium via the introduction of a temperature for 
the slow degrees of freedom that can be measured both via dynamics and 



2 



Figure 1. Sketch of the undercooled hquids motion, according to Goldstein [2]: vibrations 
in the basins of the minima of the potential energy are interspersed by more unfrequent 
change of basins. 

via statics. In Sec. 4 we speculate on the possibility of quasi-equilibrium, 
quench-rate dependent critical points. 

2. Free energy for a supercooled liquid 

The phase space of a system can be in general partitioned in different 
ways, giving rise to equivalent rewritings of the free energy [5]. In general, 
we will consider systems with an hamiltonian Ti and therefore disregard 
the momentum space as it can be trivially integrated out; we will therefore 
concentrate on the /A^-dimentional configuration space W , where / is the 
number of spatial degrees of freedoms of a particle and N is the number 
of particles. What makes us prefer or introduce a particular partition is its 
physical relevance ad interpretation. 

Stillinger and Weber introduced such a partition [3] in the spirit of the 
Goldstein picture of undercooled liquids dynamics (see Fig. 1). They asso- 
ciate to each configuration the configuration reached following a steepest 
descent path to minimize the potential energy. Such configuration is called 
an Inherent Structure (IS) and corresponds to a local minimum of the po- 
tential energy (Fig. 2). While to each IS we can associate its potential 
energy e/5, the converse is not true: to each energy level e/5 there will be 
associated a degeneracy Q{eis) counting the number of local minima with 
energy e/5 [6] (Fig. 2). Following the reasoning of Goldstein we can already 
argue that e/5 is a good candidate to be a slow variables of the system, as 
it is associated to the global structure of a liquid and not to its vibrational 
excitations. Defining the configurational entropy Sconf = kBln-{Q{eis)), we 
can write the Helmotz free energy of the system at temperature T and 
volume V as [3] 

F{T, V) = eis{T, V) - Ts,onf{eis, V) + /„fe(e/5, T, V) (1) 

where e/5 is the energy of the typical minima at (T, V), Sconf for each fixed 
V will depend implicitly on the temperature via the T dependence of e/5 
and fyib is the average free energy of the typical basins at (T, V). At fixed 
y, fvib will depend on the shape of the basins (we assume that the shape 



Figure 2. Sketch of the StiUinger- Weber partition of the configuration space W. The 
configurations indicated as /5" are local minima of the potential energy; all the configu- 
rations of the basin W axe associated to IS* via steepest descent. 

of the basin can be parametrized by ejs) and on the average kinetic energy 
(therefore on T). The vahdity of such approximations has been carefully 
checked and confirmed via numerical simulations of different atomic and 
molecular model systems [7-12] 

For real systems Sconf has been often approximated by the difference of 
the entropy of the liquid and the entropy of some reference solid phase [13]. 
For computer liquids we have access to the IS configurations and can there- 
fore extimate the "true" Sconf- In particular, we can use thermodynamic 
integrations using the relation 

deis , dfmb deis A /„s 

for the differential of Sconf obtained using the extremum condition 6F = 0. 
To calculate F, first an accurate estimate of ejs must be obtained by mini- 
mizing the potential energy of equilibrium configurations of the liquids. The 
estimate of f^ib includes two steps: first, the harmonic contribution Z'* valid 
at low temperatures is obtained calculating the normal mode frequencies 
(i = l-.fN [14]) near the IS minima; the exact calculation for a collection 
of harmonic oscillators gives us the estimate [15] 

where the average << .. >> [16] is intended both on the IS configurations 
and on all the fN modes. For systems like binary mixture Lennard Jones, 
is already a good estimate for f^n, [7]. In the case of a molecular liquid, 
the SPC/E model for water [17], measurable deviations are found [8]. Still, 
deviations are small and can be estimated as perturbative corrections to 
the ideal harmonic oscillator behavior [18]. In practice, what it is found is 
that heating up IS's one Ccin estima-te the a,nh.ctrmonic corrections Ganharm — 
aT^ + bT^ + ... to the ideal harmonic behavior e = ejs + (//2)A:^r. 
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Figure 3. The free energy of undercooled liquids can be expressed in terms of the energy 
levels eis of the inherent structures IS, to their degeneracy Q{eis) and to the shape (in 
harmonic approximation, to the curvature cu) of the associated basins. 

Therefore, once eis and fmb have been estimated, we can recontruct 
Sconf by thermodynamic integrations. The behavior of Scon/ is related to 
the possibihty of an ideal glass transition [13]. If Sconf becomes zero at 
some finite temperature Tx (the so-called Kauzmann temperature [19]), 
this means that the liquid is in a metastable state, trapped in the basin 
of a single minima; therefore diffusivity is zero and the thermal vibrations 
around the disordered structure (the IS of the minimum) are the only sort 
of motion taking place. A precise functional dependence has been proposed 
by Adam and Gibbs [20] among Sconf (a static quantity) and the diffucion 
coefficient D (a dynamic quantity): 

A 

D oce ^"-°"/ (4) 

where A is almost constant with T [20, 21]. Up to now, the Adam-Gibbs 
relation (Eq. 4) has always been found to be consistent with simulation 
results [8, 9, 11, 12] (see Fig. 2). 

3. Free energy for quasi-equilibrium states 

Attempts to construct thermodynamics for out-of-equilibrium states have 
not yet come to a success; in particular, for glassy systems the history of 
the system cannot in general be disregarded so that no static theory with 
a small number of parameters can be constructed [5]. In the case of un- 
dercooled systems, where a separation between two different time scales 
becomes sharper and sharper, it is possible to generalize in a straight- 
forward way the expression for the free energy. In the case of simulated 
liquids, quenching at low T, the time evolution of ejs is slow, with a decay 
compatible with a logarithmic or a very small power-law decay [7], while 
the short time dynamics has a much faster timescale and is vibrational in 
character (see Fig. 3). We can then treat our system as two subsystems in 
contact: a vibrational subsystem, linked to the fast vibrations in a basin, 
and the configurational subsytem of the local minima /S"s, whose dynam- 
ics is reflected by the slow decay of e/5 (see Fig. 3). On the timescales 
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Figure 4- Adam-Gibbs plot of difFusivity D versus 1/Tsconf [8]. The plot shows that 
there is a strong correspondence among the dynamics and the statics of undercooled 
liquids [22, 23]. For the sake of clarity, the values of D have been shifted by multiplying 
by powers of 5 from 5''(the lowermost curve p — 0.95g/cm'^) to 5® (the uppermost curve 
p = lAOg/cm^) 

where the vibrational subsystem comes in equilibrium with the external 
heath-bath at temperature Tf,ath (in this sense, the vibrational subsystem 
is canonical [24]), ejs is almost constant (in this sense, the configurational 
subsystem is almost microcanonical [24]). 

If we assume that the basins corresponding to a given ejs visited in 
out-of equilibrium are the same (or have the same characteristics) of the 
ones visited in equilibrium, we can use the expression of fvib{^is,T) ob- 
tained in equilibrium to take account of the vibrational part of the free 
energy. Obviously, the vibrational temperature will be equall to the exter- 
nal bath temperature Tf^ath- The configurational part corresponds to the 
almost microcanonical subsystem of the energy levels ejs together with 
their degeneracies Q{eis)- Now, we can write in general the free energy of 
this subsystem as eis — Tint s con fi^ is) where Tint (we will indicate Tint as 
the internal temperature) can be in general different from Tf^ath [5]: 



Tint interpreted as a parameter that "optimizes" the free energy [26] 
given the constraints that the energy of the minima is ejs and that the 
vibrational subsystem has temperature Ti,ath- Therefore, we can derive an 
expression of Tint [7, 27] extremizing the generalized free energy: 



F = CIS - Tints con f{e is) + fvibi^IS, Thath)- 



(5) 



deis 



dF 



Tint {e IS, Thath) 




(6) 
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Figure 5. (a) The typical decay of a correlation function in undercooled liquids. The fast 
relaxation corresponds to the vibrational dynamics and can be treated in the harmonic 
approximation using the istantaneous normal mode approach (b) Typical decay of eis 
after a deep quench. A slow logarithmic law is compatible with numerical data [7, 25]. 
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Figure 6. Quenching the system (i.e. lowering rapidly the external temperature Tbath 
of the heath bath), vibrational degrees of freedom come rapidly at a canonical (NVT) 
equilibrium with the bath, while the subsystem corresponding to structural changes (the 
/S"s) relaxes slowly and can be treated as an almost microcanonical (NVE) system in 
which eis does not change over the characteristic vibrational timescales. 

Up to now, we have a theoretical prescription for Tint that is self con- 
sistent with our writing of the generalized free energy. At this point, two 
questions naturally arise: first, is it possible to write the free energy in terms 
of experimentally measurable quantities? And then, can we measure Tj„t? 

Regarding the first problem, we recognize that already in the 30's an 
approach to the construction of an out-of equilibrium free energy for ex- 
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Figure 7. Fluctuation-dissipation plot of the response < A>p oi the system to a small 
perturbation versus the equilibrium, unperturbed decay of the cross correlations Cab 
among observables A and B. According to Eq. 7, the slope of the plot gives a measure 
of the temperature of the system. While at short times (i.e. high values of Cab) the 
temperature of the system is in equilibrium with the bath temperature, at long times 
(high values of the response < A >p) the system reveals that slow degrees of freedom 
are at a temperature higher (lower slope of the plot) than the bath. 



perimental glasses had been attempted by Davies and Jones [28] in terms 
of "frozen observables" , i.e. quantities that mantain values typical of high 
temperatures after quenching; a successful implementation of such ideas is 
found in [29] , who uses the "more" concrete entalpies of shock compressed 
hard-sphere liquids [30] to calculate the equilibrium configurational entropy. 

Regarding the measurement of Tint , we can resorts to dynamics measur- 
aments: it is in fact known both from theory [31], numerical simulations [7] 
and experiments [32], that out of equilibrium systems can exhibit at least 
two temperatures. A simple way to measure such temperatures is through 
the application of linear response theory [33]. Let's consider a pertubed 
hamiltonian Hp = H — XB, where H is the unperturbed Hamiltonian, B 
is a generic observable and A is a small parameter. Linear response theory 
predicts [34] that, switching on the perturbation at time t = 0, we have 

< A{t) >p= ^ [CABit) - CABm (7) 

where < A(t) >p (the response of the system) is the time evolution of the 
average of an observable A in presence of the perturbation and CAsit) =< 
A{t)B{0) > is the correlation function (the decay of fluctuations) between 
the observables A and B in absence of the perturbation. Therefore, from 
the slope of a parametric plot of < ^ >p versus XCabi the temperature is 
measured from the dynamics of the system (see Fig. 3). 
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While at short times such temperature coincides with Tbath, at long 
times a temperature Tg// 7^ Tf,ath is measured in quenching experiments. In 
the case of binary-mixture Lennard Jones system, it has been verified that 
Teff = Tint{eis,Tbath) [7]- In general, Tg// tends to be higher than Tbath 
and of the order of the temperature at which the system is in equilibrium 
exploring basins of depth ejs- Somehow, the configurational part of the 
system "remembers" the temperature at which was originated. 

4. Out of equilibrium phase diagram for SPC/E water 

Equation 5 and 6 allow to develop, using the inherent structure formalism, 
an expression for the free energy for disordered materials even in out-of- 
equilibrium states; for such free energy, the "history" of the system is de- 
scribed by its ejs value. In this section we present a preliminary application 
to the case of SPC/E water. 

The low temperature of SPC/E water has been object of a massive com- 
putational investigations in order to investigate the slow dynamics of molec- 
ular liquids [35] and to check the liquid- liquid critical point scenario [36] . In 
particular, the properties of the energy landscape have been throughouly 
analized in order to find connections between the dynamics and the statics 
of such a system [8, 22, 37]. 

The profile of the energy of SPC/E versus the volume develops a change 
in concavity already in the region numerically accessible (see Fig. 4), sug- 
gesting the possibility of a phase transition hindered by the entropic terms. 
Low temperature extrapolations predict indeed a liquid-liquid critical point 
a temperature of ~ 130 K and a specific density of ~ 1.1 gjcrr? [38]. Such a 
transition would be located beyond the predicted Kauzmann line for SPCE 
water [8] and therefore unaccessible to equilibrium studies [39]. 

To derive an expression for the free energy in the low temperature region 
of SPC/E water, we select for simplicity the harmonic approximation, where 

UiMs,V,T^,^th) = -kBTbath « In » (Eq. 3). We use 

the data for Sconf at T = 210 [8] as reference entropy to obtain an expression 

for Sconf- 

,h . _ ,T=2iO(y. , /^(^^^'^^ deis + df^ 

Sconf{(iIS,V) - Sconf ) + 7f , (oj 

J210 

where T^ejs, V) is the temperature at which a system at volume V in 
equilibrium has 75"s of energy ejs- We obtain then the generic expression 
for the free energy in harmonic approximation: 

F'^ieis, Tbath, V) = eis-Tint{eis, T^ath, V)s^onfieiS, ^)+/''(e/5, Thath, V). 

(9) 
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Figure 8. ejs, Sconf and fvib contributions to the free energy of SPC/E water at 
T — 210 K. The shape of e/s indicates that the system would phase separate if it was 
not for the strong —Tsconf term which avoids the change of concavity in the free energy. 
For SPC/E water, fyn, is pratically irrelevant to change the concavity of the free energy. 



First, we check the equihbrium phase diagram for F . No phase tran- 
sition is found in the physical region where Sconf > 0. We then perform an 
ideal experiment: we consider the free energy of state points in equilibrium 
at a temperature Tgq after a sudden quench at T^ath = WOK (well below the 
Kauzmann locus). After a rapid quench, ejs has not time to change on the 
timescales of vibrations, so that eis and Sconf remain to their equilibrium 
values, while /,„i5 « f^{eis{Teq), T^ath = WOK). 

The value of the internal Tint can be calculated from Eq. 6 and will 
be in general different from the temperature Tgq at which the quench has 
started. In the harmonic approximation, Eq. 6 becomes 

T^nt = Teq^-=^^ (10) 

where a = ^«^Bln^{eis,V)\» j,^^ ^-^^ temperature at which minima 
of depth eis are populated in equilibrium. 

In the case of SPC/E water, as can be inferred from Fig. 4, a is small 
(order W~^K~^) such that Tint ~ T^q, differing from the case of the binary 
mixture Lennard-Jones [7] where a is relevant and Tint < Tgq significantly. 

Therefore, for SPC/E water in the harmonic approximation, it is not 
possible to "shift up" the critical point above the Kauzmann line. This is 
true for two reasons: first, in SPC/E water the vibrational part fyn, of the 
free energy has a very small depenence both on V and on ejs compared to 
the remaining term. Therefore, if we ignore fyn,, the phase diagram stays 
almost unchanged. The only way then to get a phase transition is to "kill" 
the entropic term TintSconf'i but from Eq. 6 we see that, as dfyib/dejs « 1, 
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Figure 9. (a) example of the linear dependence of eis versus l/T (b) the average 
logarithm of the basin's curvature is approximatively linear in eis 



Tint ~ Teq where Tf.q will be approximatively equal to the temperature at 
which the quench is started as the change of eis is logarithmically slow. 
So, for systems where fmb does not vary much with eis (i-e- the shape of 
the basins does not vary much), the temperature of the slow degrees of 
freedom of the systems after a quench stays almost the same as the start- 
ing temperature of the quench. In order to decrease drastically the internal 
temperature (and therefore decrease the entropic term of the free energy 
hindering eventual critical points), one has to resort is a system with a 
strong dependence of the shape of the basins with eis- In the extreme case 
dfvib/deis >> 1, we would have the ideal situation Tint ~ T^ath- The pos- 
sibility of the appearence of critical points in out-of-equilibrium situation 
is a new issue that deserves to be further addressed and investigated. 
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